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Abstract 

Background: Freezing provol<es severe yield losses to different fall-sown annual legumes. Understanding the molecular 
bases of freezing tolerance is of great interest for breeding programs. Medicago truncatula Gaertn. is an annual temperate 
forage legume that has been chosen as a model species for agronomically and economically important legume crops. 
The present study aimed to identify positional candidate genes for a major freezing tolerance quantitative trait locus that 
was previously mapped to M. truncatula chromosome 6 (Mt-FTQTL6) using the LR3 population derived from a cross 
between the freezing-tolerant accession F83005-5 and the freezing-sensitive accession DZA045-5. 

Results: The confidence interval of Mt-FTQTL6 was narrowed down to the region comprised between markers 
MTIC153 and NT6054 using recombinant F7 and Fg lines. A bacterial-artificial chromosome (BAG) clone contig map 
was constructed in an attempt to close the residual assembly gap existing therein. Twenty positional candidate 
genes including twelve G-repeat binding factor (GBF)/dehydration-responsive element binding factor 1 (DREBl) 
genes were identified from BAG-derived sequences and whole-genome shotgun sequences (WGS). CBF/DREBl 
genes are organized in a tandem array within an approximately 296-Kb region. Eleven CBF/DREBl genes were 
isolated and sequenced from F83005-5 and DZA045-5 which revealed high polymorphism among these accessions. 
Unique features characterizing CBF/DREBl genes from M. truncatula, such as alternative splicing and large tandem 
duplication, are elucidated for the first time. 

Conclusions: Overall, twenty genes were identified as potential candidates to explain Mt-FTQTL6 effect. Their future 
functional characterization will uncover the gene(s) involved in freezing tolerance difference observed between 
F83005-5 and DZA045-5. Knowledge transfer for breeding improvement of crop legumes is expected. Furthermore, 
CBF/DREBl related data will certainly have a large impact on research studies targeting this group of transcriptional 
activators in M. truncatula and other legume species. 

Keywords: Freezing tolerance, Medicago truncatula, Quantitative trait locus (QTL), Fine mapping, Gandidate genes, 
CBF/DREBl genes. Tandem duplication. Sequence polymorphism 
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Background 

Freezing is an adverse abiotic factor with severe negative 
impacts on plant health and productivity [1]. Most tem- 
perate plants can increase their freezing tolerance after 
exposure to low temperatures, a process called cold ac- 
climation [1-3]. CBF/DREBl genes play a key role in the 
regulation of the transcriptome during cold acclimation 
[4]. They belong to the AP2/EREBP family of trans- 
cription factors and were first isolated from Ambidopsis 
thaliana (L.) Heynh. [5-7]. Up to now, CBFIDREBl 
genes have been identified in numerous herbaceous and 
woody plant species [8-11] and different studies have re- 
ported their significant role in freezing tolerance [12-16]. 

Grain legumes are an important source of food and 
feed worldwide. With their seeds containing 20% to 30% 
protein, they largely contribute to the protein require- 
ments of humans and animals [17]. Forage legumes are 
a valuable source of protein, fiber and energy for live- 
stock as well [17]. Besides their nutritional importance, 
legumes are an exceptional component to sustainable 
agriculture. In fact, most legume species are able to es- 
tablish nitrogen-fixing symbioses with rhizobial bacteria 
[17,18] which reduces the need for fertilizer application 
in crop rotation systems. Nowadays, the susceptibility to 
low temperatures and freezing conditions still represents 
a major constraint to the cultivation of legumes in many 
agro-ecological zones. The identification of the genes re- 
sponsible for the natural freezing tolerance variation can 
enhance the breeding progress and subsequently the re- 
lease of novel freezing-tolerant cultivars. 

Medicago tmncatula is an annual legume of Mediterranean 
origin [19]. It is not only a valuable forage crop as in 
Southern Australia [20] but has been chosen as a model 
molecular-genetic system for legume biology [21]. A large 
array of genomic tools has been developed for the 
reference accession A17 [22] and a draft sequence of the 
euchromatic portion of its genome comprising approxi- 
mately 94% of all genes has been released [23]. These re- 
sources can accelerate the identification of the molecular 
determinants of various traits in M. truncatula including 
freezing tolerance. Knowledge can be then transferred to 
different agronomically and economically important crops 
considering their phylogenetic closeness to M. truncatula 
and the remarkable synteny between their genomes and 
that of M. truncatula [24-28]. 

Studies on cold acclimation and freezing tolerance in 
M. truncatula are scarce [29-34]. Few biochemical changes 
have been reported in response to low non-freezing 
temperatures for this species [30,31,34]. At the molecular 
level, transcripts corresponding to M. truncatula cold 
acclimation-specific gene 15, a candidate CBFIDREBl tar- 
get, have been found to accumulate 6 hours after expo- 
sure to low temperature [32]. Four CBFIDREBl genes, 
i.e. MtCBFl-4, have been reported to be rapidly induced in 



response to cold stress [32,34,35] . Furthermore, the overex- 
pression of MtCBF3 (also named MtDREBlC) has been 
shown to improve the freezing tolerance of transgenic 
M. truncatula accession Jemalong [14]. MtCBFl and 
MtCBF4 are located on M. truncatula chromosomes 5 and 
1, respectively, but no information on the physical 
localization of MtCBF2-3 has yet been obtained. Apart 
from MtCBFl -4, no other CBFIDREBl genes have been 
described and characterized from M. truncatula even 
though expressed sequence tag (EST) and complementary 
DNA (cDNA) potentially representing transcripts from 
other putative CBFIDREBl genes are available [36,37]. 

Recently, Avia et al. [38] have uncovered natural vari- 
ation for freezing tolerance among 15 M. truncatula ac- 
cessions after a cold acclimation period. The genetic bases 
for freezing tolerance variation between the contrasted ac- 
cessions F83005-5 and DZA045-5 have been identified 
using a quantitative trait locus (QTL) mapping approach 
[38]. A major freezing tolerance QTL (Mt-FTQTL6) ac- 
counting for 40% of the phenotypic variation has been 
mapped to a region of M. truncatula chromosome 6 [38] 
coinciding with an assembly gap in the A17 euchromatic 
sequence [39]. The present study aimed to fine map 
Mt-FTQTL6 and to identify genes located in the corre- 
sponding genomic region. 

Results 

Freezing tolerance evaluation of recombinant lines 
validates Mt-FTQTL6 position and narrows down its 
confidence Interval 

Eighteen Fy or Fg lines were evaluated for freezing toler- 
ance after a cold acclimation period (Figure 1). Based on 
markers used in this study (Additional file 1), screened 
lines correspond to 9 different haplotypes; 7 of which 
carry recombination events within or next to Mt-FTQTL6 
confidence interval (Figure 1). Phenotypic data from 
F83005-5 and DZA045-5, included as check plants, were 
significantly different as expected. Comparisons of the 
phenotypes and marker genotypes validated Mt-FTQTL6 
confidence interval previously located between markers 
NT6001 and NT6019 [39]. Freezing damage scores from 
the self-pollinated progenies of lines 76-08-04 and 76- 
06-187-02 showed that F83005-5 inserts upstream 
NT6001 have no phenotypic effect. The same conclusion 
was reached concerning F83005-5 inserts downstream 
NT6019 based on the progenies of line 76-05-06. Further- 
more, Mt-FTQTL6 interval could be delineated to a 
smaller region of 0.4 cM between markers MTIC153 and 
NT6054 (Figure 1 and 2; Additional file 2). A critical re- 
combination event in line 76-02-86, as seen in lines 76-02- 
86-07 and 76-02-86-12, suggested MTIC153 as the left 
marker for Mt-FTQTL6 confidence interval. The right 
marker, NT6054, was identified based on information 
from the progenies of lines 76-09-04 and 76-11-108-04. 
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Figure 1 Phenotypic evaluation of different haplotypes carrying or not recombination events within/next to tfie confidence interval of 
Mt-FTQTL6. F7 and Fg seedlings were obtained from the self-pollination of the homozygous F5 and F7 lines whose code names are indicated in 
the 'Mother plant identifier' column. These lines originally derive from a cross between M. truncawla accessions F83005-5 [freezing-tolerant] and 
DZA04S-5 [freezing-sensitive] (see 'Development of plant material' in Methods). The genotypes of Fy/Fg lines at Mt-FTQTL6 are provided based on 
information from 26 markers that are shown according to their established genetic and/or physical order ([39]; this study). Green bars refer to 
chromosomal regions harboring alleles from F83005-5. Orange dotted bars represent regions with alleles from DZA045-5. Markers bordering the 
confidence interval of Mt-FTQTL6 according to [39] are underlined. These are 5.7 cM distant [39]. Haplotype groups A-l contain the recombinant 
and non-recombinant F/Fg lines that were compared between each other. The mean freezing damage score (FDS) of each line is provided in 
the 'Freezing tolerance assessment' part of this figure. Pairwise comparisons of mean FDS were conducted between recombinant lines, respective 
controls (if any), F83005-S and DZA045-5 (see 'Evaluation of freezing tolerance' in Methods). Letters (a) and (b) are used to distinguish significantly 
different FDS in each haplotype group. Due to the lack of individuals representing the parental accessions, the mean FDS of the progeny of 
recombinant line 76-1 1-108-04 in group I was only compared with those of the progenies of 76-1 1-108-07 and 76-1 1-51-08, considered as 
positive and negative controls respectively. Deduced phenotypes of Fy/Fg lines are indicated to the right of the figure: 'S', sensitive to freezing; 
T, tolerant to freezing. Grey-highlighted markers in the 'Genotype at Mt-FTQTL5' part are those contained in Mt-FTQTL6's confidence interval 
delimited according to genotype and phenotype data from screened lines. 



Full-BAC sequencing and WGS assembly provide 
sequences spanning the assembly gap coinciding with 
Mt-FTQTL6 

Efforts have been made to close the assembly gap between 
markers NT6001 and NT6019. Six primary BAC clone 



contigs partly spanning Mt-FTQTL6 could thus be con- 
structed (Additional files 3, 4, 5). Mt-FTQTL6 confidence 
interval delimited through fine mapping (see above) is 
partly covered by BAC clone contig IV (Additional file 5). 
In order to generate the candidate genomic sequence of 
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Figure 2 Genetic and physical map positions of Mt-FTQTL6 candidate genes on M. truncatula chromosome 6. Genetic markers shown at 
the top of the figure refer to gene-based markers bordering or located in the confidence interval of Mt-FTQTL6 (refer to Figure 1 and Additional 
file 2). BAC clones originating from the genomic region corresponding to Mt-FTQTL6 are presented according to their relative positions and with 
respect to genetic markers (see Additional file 5 for details on the BAC clone contig map). The identifiers of the BAC clones whose inserts were 
sequenced in this study are underlined. Scaffold065bis is a WGS scaffold from A17 that largely covers the confidence interval of Mt-FTQTL6. 
Arrowheads indicate the position and the transcriptional orientation of candidate genes. No information on the transcriptional orientation of 
MtBAG-1 is available, a, MtBAG-1; b, MtCBFM; c, MtCBF13; d, MtCBFS; e, MtCBFU; f, MtCBFl 1; g, MtCBF2; h, MtCBFW; i, MtCBF9; j, MtCBFS; k, MtCBF7; 
I, MtCBF6; m, MtCBFS; n, MtPERLD; o, MTR_050s0019; p, MTR_050s002a, q, MtZFWD; r, MTRJJSAsOOOl; s, MTR_054s0019, t, MTR_6g089580. 
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Mt-FTQTL6, inserts from BAC clones mth2-92015, mtel- 
60A22 and mth2-221P20, were first fully-sequenced 
(GenBank accession numbers [GB acc] KF006382-84). Se- 
quence information from these 3 clones (Figure 2) extended 
BAC clone mth2-53124 insert sequence (GB acc AC229695) 
with 217,951 bp assembled in 17 contigs (Additional file 6). 
With Illumina WGS sequences becoming available from 
A17 (F. Debelle, personal communication), BAC by BAC se- 
quencing has not been kept up. Instead, marker and BAC- 
end sequences associated with Mt-FTQTL6 were used to 
search for scaffold sequences corresponding to the target 
region. A scaffold of more than 1 Mb long, namely scaf- 
fold065bis, was thus identified. It covers BAC clone con- 
tigs IV and V (Additional file 5) and extends mth2-221P20 
insert sequence with approximately 155 Kb (Figure 2). 

Twenty non-transposon genes reside in IVlt-FTQTL6 
confidence interval 

Twenty non-transposon genes with significant matches to 
GenBank's non-redundant nucleotide (nr/nt), expressed 
sequence tag and protein databases, including those corre- 
sponding to markers MTIC153 and NT6054, could be 
identified in the candidate sequence for Mt-FTQTL6 
(Figure 2; Table 1). Spacing between adjacent genes varies 
from 1,586 bp (between MtCBFS and MtPERLD) up to 
318,549 bp (between MTR_054s0019 and MTR_6g089580). 
Coding sequences have an overall G + C content of 
42.15%. Sixteen genes are supported by cognate EST and 
additionally 7 of them by cDNA sequences (Table 1). Com- 
parisons between genomic and EST/cDNA sequences in- 
dicate that four genes (i.e. MtCBFS, MtCBF6, MtPERLD 
and MTR_054s0019) have 2 transcript isoforms each. 

All positional candidate genes for Mt-FTQTL6 were 
found to have highly similar counterparts across a set of 
dicotyledonous species for which genome sequence infor- 
mation is available (Additional file 7). AU, except for hom- 
ologous genes to MTR_054s0001, are contained in 
syntenic blocks exhibiting large gene order conservation 
with Mt-FTQTL6 (Additional file 7). Three chromosomal 
segments in A. thaliana and in Solarium lycopersicum L., 
two segments in Populus trichocarpa (Torr. & Gray) and a 
single segment in Vitis vinifera L. share colinearity with 
Mt-FTQTL6. Among legumes, two additional M. tmncatula, 
six Glycine max (L.) Merr. and three Lotus japonicus 
L. regions are colinear with Mt-FTQTL6. Few rear- 
rangements could be noted among conserved regions 
including gene loss, local duplication or also transloca- 
tion events (Additional file 7). 

Twelve candidate genes organized in a tandem array 
belong to the CBF/DREB1 group of the AP2/EREBP 
transcription factor family 

CBFIDREBl genes represent 60% of all candidates for Mt- 
FTQTL6 (Table 1). They are organized in a tandem cluster 



spanning a region of approximately 296 Kb (Figure 2). No 
intervening spacer genes, except those for retrotranspo- 
sons, exist in between. MtCBF12 and MtCBFlS are pseu- 
dogenes. A deletion of a guanine at position 314 in 
MtCBF12 coding sequence causes a frameshift and a pre- 
mature stop codon at the 148* amino acid. MtCBFlS cod- 
ing region shows an in-frame 848-bp insertion and two 
point mutations leading each to a premature stop codon. 
MtCBF2-3;5-ll;14 are likely to represent functional 
protein-coding genes. The corresponding proteins share 
61 to 92% overall amino acid identity (Additional file 8) 
and all show typical features of the CBF/DREBl group 
(Figure 3). They have each a 58-amino acid long AP2/ERF 
nuclear signaling and DNA binding domain [42-44] and 
an acidic C-terminal region [isoelectric point (pi) ranging 
between 3.94 and 5.04] that may function in trans- 
activation [45]. MtCBFll contains in addition to the 
58-amino acid AP2/ERF domain, a 17-amino acid long 
sequence most likely corresponding to the C-terminal end 
of a second AP2/ERF domain (Figures 3 and 4). CBF signa- 
ture sequences comprising PKKP/RAGRxKFxETRHP and 
DSAWR motifs bracketing the AP2/ERF domain [8] are 
largely conserved among MtCBF2-3;5-ll;14. These latter 
also show A(A/V)xxA(A/V)xxF [46] and LWSY motifs 
[47], reported to be conserved in CBF/DREBl proteins in 
different plant species. Variants at the CBF/DREBl charac- 
teristic amino acid sequences are listed in Additional file 9. 
The (L/Y)(L/Y)x(N/S)(M/L)A(E/Q)G (M/L)(L/M)xxPP se- 
quence, previously suggested as a CBF/DREBl conserved 
motif [46] is not present in MtCBF2-3;5-14 except for the 
NMA motif. These proteins have rather a NM(A/V) 
LMSPTHS conserved sequence at the same position. 

Distinct levels of polymorphism exist between 
MtCBF2-3;5-12;14 alleles from F83005-5 and DZA045-5 

MtCBF2-3;5-12;14 specific primers yielded positive amp- 
lification from F83005-5 and DZA045-5 indicating that 
all 11 genes are present in these accessions. No amplifi- 
cation could be obtained with MtCBFlS-specific primers 
from both accessions despite using different primer 
combinations and PCR conditions. For MtCBF2, two 
distinct PCR products were obtained from F83005-5 
(named MtCBF2a and MtCBF2b) using distinct forward 
primers but not from DZA045-5, thus indicating that this 
gene has undergone duplication at least in F83005-5's 
background. Overall, 14,953 bp (F83005-5) and 15,038 bp 
(DZA045-5) of comparable sequences were generated 
(Additional file 10). They show a total of 402 single nu- 
cleotide polymorphism (SNP), 41 insertion-deletion (indel) 
and 9 simple sequence repeat (SSR) polymorphisms with 
206 SNP 8 indel and 3 SSR located in the coding 
sequences (Additional file 10). Contrary to A17 and 
F83005-5, DZA045-5 has an insertion of a 160-bp short 
interspersed nuclear element (SINE) in the coding 
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(1} Corresponding EST and cDNA were identified through BLAST searches against M. truncatuia EST and nr/nt databases in GenBanl< [36] (January 2013}, respectively. 1 


n case of alternate splicing, the sequence(s) 


corresponding to the less represented transcript isoform is(are) in italics; (2) similar to [41], the C-terminal domains of MtCBF2-3;5-1 1;14 


were considered as the amino 


acid sequences that occur after the signature 


motif "DSAWR" and that up to the last amino acid of the protein; (a) a cDNA sequence corresponding to MtCBFS 


was directly submitted to GenBank [36] by Chen JR, Guo L, Wang H in 2007. The respective locus was 



referred to as DREB1A; (b) In addition to MTR_050s0019 and MTR_050s0020, 7 other non-transposon genes could be identified on AC229695 but were not considered in this study. In fact, the corresponding genomic re- 
gion for these latter, spanning 114,402 kb (1-19,425; 23,144-90,266; 147,607-159,557; 185,484-201,386 bp) is completely missing from Scaffold065bis. Because of its size exceeding 200 Kb, AC229695 seems to have re- 
sulted from sequencing inserts of at least two BAC clones, one of which did not originate from Mt-FTQTL6-containing region; (c) The intron/exon structure of MTR_050s0020 from M. truncatuia Mt3.5 genome assembly 
has been revised and Its corresponding coding sequence found to contain 15 introns instead of 16. M. truncatuia EST and close protein sequences together with EST from other Papillonoideae species and notably a 
Pisum sativum L. sequence generated by RNAseq (Alves Carvalho et al. in preparation) were used for this purpose; (d) As related EST do not fully cover MtZFWD coding sequence, intron positions were predicted after 
multiple alignment with the coding sequences of Glymal6g32370 and Glyma09g27300 and EST sequences from Lathyrus odoratus L. (GO316201, GO317026, G0319782), L. japonicus (FS347722, FS346164, FS360001) 
and Phaseolus acutlfoHus A. Gray (H0781358); (e) As the EST BQ138417 does not cover the entire coding sequence of MTR_054s0001 , highly similar EST sequences from L japonicus (GO014516), Vigna unguiculata (L.) 
Walp. (FF549909} and G. max (AW1 84979, FK01 1 826) were used for intron location. The positions and sizes of introns 9 and 1 0 are different from Mt3.5. 
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Figure 3 Alignment of predicted MtCBF2-3;5-14 sequences. Hypothetical protein sequences are used for MtCBF12-13. These were obtained after 
eliminating truncating mutations from the corresponding coding sequences. For MtCBFS and MtCBFS, proteins encoded by both transcript variants are 
shown. MtCBFl 1 protein sequence, likely resulting from a chimeric gene, is separated from others. Conserved amino acid residues are in white. They 
are shaded in black when shared by all sequences or in grey if a single exception exists. Conserved amino acid residues in the C-terminal regions of 
MtCBF5;6.2;7;8.2 resulting from intron-containing transcripts are in black and grey-shaded. The AP2/ERF domain is underlined and conserved amino 
acid residues therein are in red. Open squares show conserved valine and glutamic acid residues at the 14th and 19th position of the AP2/ERF domain 
similarly to other CBF/DREB1 sequences; V14 is critical for determining the DNA binding specificity of CBF/DREB1 proteins [48]. Solid squares indicate 
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hydrophobic residues [45] that may contribute positively (full lines) or negatively (broken lines) to the trans-activating properties of MtCBF2-3;5-1 0; 
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Figure 4 Polymorphism between F83005-5 and DZA045-5 in the coding sequences of MtCBF2-3;5-12;14. The order of sequences reflects 
that of the respective genes in Mt-FTQTL6 region. Schematic representations of the coding sequences of MtCBF2-3;5-12;14 are based on the 
shortest allele (F83005-5 or DZA045-5) at each time. Blue boxes illustrate regions encoding AP2/ERF domains while orange boxes depict coding 
regions for characteristic CBF/DREB1 signature sequences, i.e. (from left to right) PKKP/RAGRxKFxETRHP, DSAWR, A(A/V)xxA(A/\/)xxF and LWSY. 
Both alternatively spliced isoforms of MtCBF6 and MtCBFS are shown. MtCBF2a and MtCBF2b, obtained from F83005-5 using primer sets designed 
on the same gene sequence from A17, are both compared to the unique sequence obtained from DZA045-5 with these primers. Black and red 
horizontal lines indicate the positions of synonymous and non-synonymous SNP, respectively. Dotted horizontal lines indicate the positions of 
indel or SSR polymorphisms. Arrows are used to delimit comparable sequences between F83005-5 and DZA045-5 in case only a partial sequence 
could be obtained from one or both accessions. Asterisks followed by letters "F" or "D" show the positions of stop codons for F83005-5 and 
DZA045-5, respectively. 



sequence of MtCBFM and another SINE element in the 
5 ' upstream region of MtCBFS. F83005-5 has a specific in- 
sertion of a putative miniature inverted-repeat transpos- 
able element in the 5 ' upstream region of MtCBF6. 

Figure 4 illustrates allelic variation between F83005-5 
and DZA045-5 in the coding regions of MtCBF2-3; 
5-12;14. As it could be noticed, nucleotide changes leading 
to differences at the amino acid level are not evenly 



distributed amongst genes. Amino acid substitution differ- 
ence between F83005-5 and DZA045-5 ranges between 1 
(MtCBFS) and 34 (MtCBFS). MtCBFM and MtCBFS are 
likely to encode truncated proteins in freezing-sensitive 
accession DZA045-5 but not in freezing-tolerant acces- 
sions F83005-5 and A17. In contrast, MtCBFU is likely to 
encode a full-length protein in DZA045-5 which is not the 
case for F83005-5 and A17. 
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Discussion 

M. truncatula is a valuable forage crop and also a prom- 
inent model for legume genomics. Freezing tolerance 
QTL have been previously identified in M. truncatula 
[38]. In the current study, Mt-FTQTL6 has been finely 
mapped by genetically and physiologically characterizing 
recombinant haplotypes carrying recombination events 
within the QTL confidence interval. Genomic data have 
been developed to fill the genome assembly gap in 
Mt-FTQTL6 region and have permitted to reveal twenty 
positional candidate genes. Interestingly, twelve genes 
belong to the CBF/DREBl group of the AP2/EREBP 
transcription factor family and occur in tandem array. 

Current knowledge on the functions of the twenty 
positional candidate genes for IVlt-FTQTL6 

The fine mapping step has significantly reduced Mt- 
FTQTL6's confidence interval and has subsequently 
allowed the identification of 20 candidate genes. Functional 
data regarding these candidate genes are still limited at 
present (Additional file 11). Only homologs/orthologs of 
MtCBF2-3;5-14, MtBAG-1, MTR_054s0019 and MtZDP 
have been studied in plant species [13,15,35,51-54]. For 
MtPERLD and MTR_050s0020 which are members of 
evolutionary conserved gene families, only non-plant 
homologs have been described [55-57]. MtZFWD contains 
both a CCCH zinc finger domain and WD40 repeats and 
thus belongs to a plant-specific subgroup of CCCH zinc 
finger protein family [58] that is not yet functionally 
described. No data are available for MTR_050s0019 and 
MTR_6g089580 to infer putative functions. Out of all can- 
didate genes (Table 1), only MtCBF2 and MtCBFS have 
been functionally characterized in M. truncatula. MtCBF2 
and MtCBFS were reported to display both a rapid and 
transient accumulation of transcripts in the leaves of 
3-week old M. truncatula plants after exposure to low 
temperature treatment [32,34]. Furthermore, transgenic 
M. truncatula lines overexpressing MtCBFS were shown 
to exhibit an improved freezing tolerance compared to the 
wild-type [14]. However, even if the CBF/DREBl genes 
can be considered as the most likely candidates for the 
Mt-FTQTL6 effect, none of the non-functionally charac- 
terized genes can be excluded at present. Further know- 
ledge needs to be gained regarding the different candidate 
genes and the QTL of interest itself. 

Co-location between a freezing tolerance QTL 
and CBF/DREBl genes in a legume species 

CBF/DREBl genes were found to co-locate with QTL for 
freezing tolerance in several plant species including 
A. thaliana [59,60], Triticum monococcum L. [61], Triticum 
aestivum L. [62] and Hordeum vulgare L. subsp. vulgare 
[63]. A large deletion in the promoter region of ^. t/ialiana 
CBF2 [59] and co-locations of freezing tolerance QTL in 



T. monococcum [61] and H. vulgare subsp. vulgare [63] with 
QTL for expression of cold-induced CBF/DREBl target 
gene C0R14b and accumulation of COR14b protein sug- 
gested CBF/DREBl genes being responsible for the QTL 
effects. The identification of CBF/DREBl genes among pos- 
itional candidates of Mt-FTQTL6 provides an additional 
example of a co-location between a freezing tolerance QTL 
and this group of genes but for the first time in a legume 
species. It will be of great interest to inspect such co- 
location in other legumes especially P. sativum having a 
freezing damage QTL syntenic to Mt-FTQTL6 [39] . 

Large size of the CBF/DREBl group in M. truncatula 

This study has permitted the physical positioning of 
twelve M. truncatula CBF/DREBl genes through the as- 
sembly of Mt-FTQTL6 region. Apart from MtCBF2 and 
MtCBFS that were identified by Pennycooke et al. [32] 
and DREBIA whose cDNA sequence was directly sub- 
mitted to GenBank (GB acc DQ778006), the other CBF/ 
DREBl genes associated with Mt-FTQTL6 were not previ- 
ously annotated or described (Table 1). As the name 
"MtCBF4" has been recently assigned to a CBF/DREBl 
gene located on M. truncatula chromosome 1 [35], the 
novel genes including DREBIA were here given consecu- 
tive numbers from 5 to 14 starting from the most prox- 
imal to MtPERLD (Figure 2). 

MtCBF2-S;5-14 occur in tandem array and are in direct 
orientation (Figure 2) which suggests that they are most 
likely derived from an unequal crossing over mechanism 
[64,65]. However, the contribution of transposable ele- 
ments to this duplication event cannot be totally excluded 
[65], especially with several transposon and retrotrans- 
poson insertions being evident between CBF/DREBl genes 
(data not shown). 

According to Young et al. [23], M. truncatula has ex- 
perienced high rates of local gene duplication compared 
to other plant genomes. It may thus not be surprising to 
have such a large tandem array organization for the 
MtCBF2-S;S-14 locus. Similar tandem duplications have 
been described for nucleotide-binding site-leucine-rich repeat 
genes [66] and lipoxygenase genes [67] on M. truncatula 
chromosomes 6 and 8, respectively. However, the main 
particularity of the tandem cluster in Mt-FTQTL6 region 
is that it concerns CBF/DREBl genes. Large groups of 
CBF/DREBl genes have always been thought to be limited 
to monocotyledons. This study reports for the first time 
that the CBF/DREBl genes can also be highly duplicated 
in legume species. M. truncatula has at least 17 
CBF/DREBl genes (considering MtCBF2a and MtCBF2b 
from F83005-5) all located in homoeologous regions 
(Additional file 7). A recent report has similarly shown 
that Eucalyptus grandis genome contains 17 CBF/DREBl 
genes [68] . It will be next important to determine if other 
legumes have as many genes that belong to the CBF/ 
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DREBl group as M. truncatula. First evidences denote 
that there may not be a general trend at this level. As 
shown in Additional file 7, G. max does not have any du- 
plication of CBFIDREBl genes in Mt-FTQTL6 syntenic 
regions. In contrast, L. japonicus has 7 tandemly-arrayed 
putative CBFIDREBl genes on chromosome 4. 

Complexity of the CBF/DREB1 locus on M. truncatula 
chromosome 6 

As for other duplication modes, tandem duplication 
generates identical genes that evolve in different ways 
[65,69-75]. In the case of the locus harboring MtCBF2- 
3;S-14, some evidences of pseudogenization are apparent: 
premature stop codons leading to truncated proteins exist 
within the corresponding coding sequences of MtCBFlS 
and MtCBF12 in A17, MtCBFU in F83005-5 and MtCBFS 
and MtCBFM in DZA045-5 (Figure 4). However, most 
CBFIDREBl duplicates are likely to encode functional pro- 
teins. Up to date, information regarding the expression pro- 
files of MtCBF2-3;5-14 is still greatly lacking. M. truncatula 
gene expression atlas [76] provides only a first glimpse 
on the expression of some genes among MtCBF2-3;5-14 
under unstressed conditions (Additional file 12). It is 
thus crucial to conduct a complete survey on the ex- 
pression patterns (including under cold stress) of all 
CBFIDREBl genes on M. truncatula chromosome 6. 
This will determine whether tandemly duplicated CBFI 
DREBl genes have evolved or not different expression 
patterns and so possibly distinct functions as already re- 
ported for different CBFIDREBl genes in other plant 
species [10,41,77,78]. The characterization of MtCBFll 
deserves a particular attention as its predicted protein 
product contains a full-length and an additional 
17-amino acid long partial AP2/ERF domain (Figure 3 
and 4). MtCBFll has most likely arisen from an in- 
frame fusion of the 5 ' portion of the coding region from 
a first CBFIDREBl gene with the 3 ' portion of the cod- 
ing sequence from a second gene. 

Alternative splicing events can diverge between dupli- 
cates, including organ- and stress-specific differences, 
and result in functional variation [79]. Based on genomic 
sequence alignment with EST and cDNA data, MtCBFS 
and MtCBFS have two transcript variants each, gen- 
erated by intron retention/splicing events (Table 1). 
MtCBFS and MtCBF7 contain each an intron at the 
same position as for MtCBFS and MtCBFS intron- 
spliced isoforms but no evidence yet suggests that 
MtCBFS and MtCBF7 may have intron-retaining tran- 
script forms similar to MtCBFS and MtCBFS. As EST 
databases can be incomplete with regard to alternative 
splicing transcript variants, reverse transcription PGR or 
transcription analysis using second generation sequencing 
would allow determining alternative splicing patterns with 
more precision. This is the first study to provide data on 



the presence of introns and alternative splicing forms in 
CBFIDREBl genes. Each of the isoforms oiMtCBFS-8 will 
need to be characterized apart in order to determine if 
these are functional and then compare them to other 
CBFIDREBl genes and alternatively spliced forms if any. 
Any possible association of alternative splicing of CBFI 
DREBl genes with M. truncatula response to low temper- 
atures should be inspected. Recent reports from grasses 
have shown the existence of alternatively spliced isoforms 
for AP2/EREBP family members [80-82] homologous to 
A. thaliana DREB2 [6] which is known to be activated by 
dehydration and high-salt stresses. The alternative splicing 
of these genes was suggested to be important for reducing 
the activity of their target genes under unstressed con- 
ditions and thus preventing negative effects on plant 
growth. 

Conclusions 

The present study has permitted to narrow down the con- 
fidence interval of a major freezing tolerance QTL on 
M. truncatula chromosome 6 and to identify genes lo- 
cated therein. No similar studies related to legume freez- 
ing tolerance have to date been reported. Candidate genes 
represent important elements toward positional cloning of 
Mt-FTQTL6. Information gained from M. truncatula be- 
ing a model species should be transferred to agronomically 
and economically important crop legumes. The BAG 
clone contig map constructed for the assembly gap coin- 
ciding with Mt-FTQTL6 will be interesting in efforts to 
clone important loci mapped to the same region as Mt- 
FTQTL6 including a locus for resistance to Subterranean 
clover mottle virus [83] and a QTL for the number of sec- 
ondary branches [84]. When sequenced from F83005-5 
and DZA045-5, MtCBF2-3;S-12;14 coding sequences 
showed distinct levels of indels and non-synonymous 
mutations which suggest that they have evolved in a 
heterogenous manner. Allelic variation identified on cod- 
ing and immediate flanking sequences of MtCBF2-3;S- 
12;14 will be useful for future studies to establish the 
adaptative value of the different CBFIDREBl duplicates. 

Methods 

Development of plant material 

M. truncatula LR3 population comprises 178 Y-i-s recom- 
binant inbred lines (RIL) obtained by single seed descent 
from a cross between the freezing-tolerant accession 
F83005-5 and the freezing-sensitive accession DZA045-5 
[85]. F83005-5 derives from a natural population collected 
in France and DZA045-5 from an Algerian population 
[86]. Genotypic data from the LR3 population [38] indi- 
cate that twelve RIL (RIL4, 13, 17, 22, 76, 90, 101, 151, 
161, 165 and 241) have residual heterozygosity at Mt- 
FTQTL6 peak marker, MTIG153. Among them, RIL76 is 
particularly important because: (1) it is homozygous at the 
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confidence intervals of both other freezing tolerance QTL 
from the LR3 population mapped on chromosomes 1 and 
4 [38]; (2) it carries the sensitive parent's alleles at both 
these QTL; and (3) it shows a low overall genome hetero- 
zygosity. RIL76 sister lines were therefore chosen for the 
development of inbred lines showing recombination 
events at Mt-FTQTL6. Twelve F5 lines were grown and 
submitted to a genotypic analysis using: SSR markers 
evenly spaced on the eight M. truncatula chromosomes 
(3 to 6 markers per chromosome; 33 markers in total), the 
nearest SSR marker to the peak of each of the freezing tol- 
erance QTL on chromosomes 1 and 4, and SNP and SSR/ 
indel markers closely linked to Mt-FTQTL6 (6 markers in 
total). Information regarding Mt-FTQTLS-linked markers 
is available in Additional file 1. The complete list for the 
rest of markers is available upon request. Similarly to 
RIL76, sister lines did not show heterozygosity for any of 
the background markers and were all homozygous carry- 
ing alleles from the freezing-sensitive parent at the free- 
zing tolerance QTL on chromosomes 1 and 4. Five 
heterozygous recombinants at Mt-FTQTL6 were identi- 
fied, namely 76-02, 76-05, 76-07, 76-08 and 76-09. Taking 
advantage of their heterozygosity for all markers associ- 
ated with Mt-FTQTL6, 76-06 and 76-11 were selected for 
the construction of large segregating populations. Line 
76-02 was also used for this purpose because it is only 
homozygous for two markers at the upper border of Mt- 
FTQTL6, namely NT6005 and NT6001 (Figure 1). Fg 
populations that were obtained through self-pollination of 
76-02, 76-06 and 76-11 comprised 447, 241 and 232 
plants, respectively. Progenies were genotyped using 7 
SNP and 6 SSR/indel markers located on both sides of 
marker MTIC153 (Additional file 1). For recombinant in- 
dividuals, genotyping was further completed with 13 add- 
itional SNP markers (Additional file 1). A high-resolution 
linkage map constructed from a subset of marker data 
scored on 76-06 and 76-11 progenies is described in Tayeh 
et al. [39] (Additional file 2). F5 and Fg heterozygous re- 
combinant individuals were self-pollinated. Homozygous 
(Fg or F7) recombinants were selected from their progen- 
ies using corresponding markers. F7 and Fg plants were 
finally obtained through the self-pollination of homozy- 
gous recombinant plants and were submitted to freezing 
tolerance tests (Figure 1). 

Plant growth 

Seeds of M. truncatula were scarified using sand paper in 
order to rupture the seed coat. Scarified seeds were soaked 
in distilled water for 6 hours and spread on moistened filter 
papers in Petri dishes. After a 3-day dark storage at 6°C to 
break embryo dormancy and synchronize germination, 
Petri dishes were held at 20°C for three other days. Seed- 
lings were then transplanted in 2 L pots containing a mix 
of peat and compost or in pre-moistened 41 x 42-mm 



jiffy-7 pellets (Jiffy France S.A.R.L., Trevoux, France), 
depending on the experimental need. In case of non- 
germinated seeds, especially when they were recendy har- 
vested, 3 to 5 parts per million of 2-chloroethylphosphonic 
acid were applied to Petri dishes to break embryo dor- 
mancy and seedlings were transplanted 24 hours later. Ex- 
cept for freezing tolerance tests, plants were finally left to 
grow in the greenhouse at an average temperature of 18/ 
14°C (day/ night) and a 16-hour photoperiod. Seeds were 
extracted from ripe pods using in most cases a homemade 
extraction system based on a rubber mat and a plaster 
hawk as described by Garcia et al. [87]. Large extractions 
were performed using a seed mill. 

Marker development and genotyping 

In addition to markers from Tayeh et al. [39], 8 new 
markers associated with Mt-FTQTL6 were developed in 
this study (Additional file 1). Primers were designed from 
BAC-end sequences using the Primer3 tool [88] with 
default parameters. Maximum product length was 
limited to 500 bp. Total genomic DNA was isolated 
from fresh leaflets harvested on one-month-old plants 
using Nucleospin 96 Plant II (Macherey Nagel, Germany) 
according to manufacturer's instructions. Genotyping of 
the plant material from this study with background and 
Mt-FTQTL6-associated markers was conducted as de- 
scribed in Tayeh et al. [39]. Capillary gel electrophoresis 
was employed for SSR and indel markers. The high- 
resolution melting technique (in case of Fy/Fs lines) and 
the competitive allele specific PGR method (in case of 
447-plant 76-02-derived population) were used for SNP 
markers. 

Evaluation of freezing tolerance 

Freezing tolerance tests were performed in a climate- 
controlled walk-in 2.6 x 2.3 x 2 m growth chamber. Two 
experiments, consisting of two replications each, were car- 
ried out. Each replicate was composed of: (1) F7 or Fg 
plants carrying one (or two) recombination event(s) within 
Mt-FTQTL6-containing chromosomal region, (2) plants 
obtained from the same parental lines as the recombinant 
ones or from sister parental lines but that ideally do not 
have recombination events within Mt-FTQTL6 region 
and (3) parental accessions F83005-5 and DZA045-5 
(Figure 1). Ten individuals from each of the above-cited 
lines were evaluated per replicate. Seedlings, transplanted 
in Jiffy pellets (see plant growth section), were placed in 
10 X 10 well-containing homemade Styrodur plates to en- 
sure root insulation. Four different phases were applied 
during each freezing tolerance test with the following en- 
vironmental conditions: (1) nursery phase-18 days: 20°C/ 
14°C, 14-hour photoperiod and 250 [^mol.m' .s' photo- 
synthetically active photon flux density (PPFD) provided 
by fluorescent lights; (2) cold acclimation phase-20 days: 
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8°C/2°C, 14-hour photoperiod and 250 |imol.m'^.s'^ PPFD; 
(3) freezing phase-8 days: 4°C/-6°C, 10-hour photoperiod 
and 150 [imol.m'^.s'^ PPFD; and (4) recovery phase-16 days: 
16°C/5''C, 14-hour photoperiod and 250 |imol.m'^.s'^ PPFD. 
Continuous acquisition of chamber air temperature was 
carried out using a Campbell CRlOx system equipped with 
a Humitter 50Y integrated humidity and temperature 
transmitter (Vaisala, Vantaa, Finland). Temperatures at 
leaf- and root- levels were also continuously surveyed 
using eight and four homemade T-type thermocouples 
connected to the acquisition system, respectively. Photo- 
synthetically active radiation was checked twice per test 
phase by means of a LiCor 190SA quantum sensor (Licor 
Biosciences, Lincoln, Nebraska, USA). Irrigation was per- 
formed manually and stopped 2 days before and all along 
the freezing phase. Developmental stages were scored 
plant-by-plant six times during the nursery and cold accli- 
mation phases. At the 16'^ day of rewarming, plants were 
evaluated for freezing tolerance using a visual rating scale 
based on leaf injuries and ranging from 0 (no freezing in- 
juries) to 5 (all leaflets are yellow and dry: dead plant). 
Plants were scored 1 when less than one quarter of the 
leaflets were visually damaged, 2 for one quarter of the 
leaflets being affected, 3 when freezing damage was ob- 
served on half of the leaflets and 4 when three-quarter of 
the leaflets were concerned. Pairwise comparisons of dam- 
age scores from recombinant, non-recombinant and par- 
ental lines were performed by one-way analysis of 
variance followed by a Tukey test at P < 0.05 using R [89] 
(Figure 1). For each of these pairwise comparisons, the re- 
combinant line and its non-recombinant counterpart(s) 
were compared between each other and with the closest 
plants of F83005-5 and DZA045-5 within the growth 
chamber. 

Full BAC sequencing 

BAC DNA was isolated using Nucleospin fast purification 
Kit (Macherey-Nagel, Diiren, Germany). Full BAC sequen- 
cing was performed using 454 multiplexing technology on 
a GS-FLX Titanium sequencer (454 Life Sciences, Roche 
Diagnostics, Branford, Connecticut, USA). Raw reads were 
cleaned using PyroCleaner [90]. Reads sharing more than 
95% homology with BAC vector or Escherichia coli gen- 
ome sequence were excluded, de novo assembly of the 
cleaned reads was performed using Newbler (version 2.3; 
454 Life Sciences, Roche Diagnostics). The resulting 
contigs were, when possible, ordered and oriented based 
on marker sequences located therein. End sequences 
from overlapping BAC clones (Additional file 5) were 
used for the same purpose. Complementarily, end se- 
quences of the different contigs were submitted to 
blastn/blastx [91] search against M. truncatula nr/nt, 
high-throughput genomic sequence, EST, genome sur- 
vey sequence and protein databases in GenBank [36] 



at the National Center for Biotechnology Information. 
End sequences sharing more than 99% identity with 
the same target sequence(s) were considered as indi- 
cative of adjacent contigs. 

Sequencing of CBF/DREB1 genes from F83005-5 and 
DZA045-5 

The isolation and sequencing of MtCBF2-3;5-14 from 
F83005-5 and DZA045-5 was undertaken in order to de- 
termine if all genes composing the CBF/DREBl cluster in 
A17 are similarly present in these accessions and to iden- 
tify the polymorphism in the coding and immediate flank- 
ing sequences of these genes. Two to 3 specific forward 
and also 2 to 3 specific reverse primers were designed per 
gene. PCR were carried out using a touchdown protocol 
with the following thermal cycling conditions: 4 minutes 
at 94°C; 5 cycles at 94°C for 30 seconds, annealing 
temperature (Ta) [-l°C/cycle] for 30 seconds, 72°C for 
90 seconds; 35 cycles at 94°C for 30 seconds, [Ta-5°C] for 
30 seconds, 72°C for 90 seconds; and 10 minutes at 72°C. 
Different primer combinations and PCR conditions 
(Ta and/or MgCl2 concentration) were tested to optimize 
the amplification result. PCR products from the optimal 
conditions were purified from agarose gel using Nucleos- 
pin Extract II (Macherey-Nagel, Diiren, Germany) follow- 
ing manufacturer's instructions and Sanger-sequenced 
using BigDye Terminator v3.1 chemistry (Applied Biosys- 
tems, Foster City, California, USA) on a 3130x1 Genetic 
Analyzer (Hitachi/ Applied Biosystems). Primers and amp- 
lification conditions that were used in the sequencing step 
are provided in Additional file 10. For 9 out of the 12 
CBF/DREBl genes, one or two internal primers were 
needed to ensure full sequence coverage and/or overcome 
sequencing difficulties caused by simple sequence repeats 
[(A)n or (T)n] present in 5'- or 3 ' -non-coding flanking re- 
gions (Additional file 10). Base calling was performed with 
Sequencing Analysis Software (version 5.4; Applied Bio- 
systems). Electropherograms were manually trimmed of 
poor sequence data. 

Availability of supporting data 

The data sets supporting the results of this article are 
included within the article and its additional files. All se- 
quences were submitted to GenBank [36]. Insert sequences 
of BAC clones mth2-92015, mtel-60A22 and mth2- 
221P20 were deposited in the "high throughput genomic 
sequences" division under accession numbers KF006382- 
84. BAC-end sequences were submitted to the "genome 
survey sequences" division and have been assigned the 
accession numbers JY974377-472. CBF/DREBl gene 
sequences from A17, F83005-5 and DZA045-5 were 
deposited in the "plant, fungal, and algal sequences" 
division with accession numbers KC997199-225. 
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Additional files 



Additional file 1: List of markers used for plant material 
development and BAC clone contig map construction. Provides a list 
of the markers that were used for one or more of the following 
purposes: (1) genotyping of segregating populations; (2) genotyping 
of F7/F8 lines used for Mt-FTQTL6 fine mapping; (3) synthesis of the PCR 
products used as probes for BAC filter screening and (4) validation of 
positive BAC clones and BAC contig construction. It includes primer 
sequences, PCR conditions and information on the type of polymorphism 
for corresponding markers. 

Additional file 2: Genetic linkage map of Mt-FTQTL6 region. 

Contains a partial linkage map of M. truncatula chromosome 6 which 
was constructed using data collected on 453 Fs plants. It shows the 
relative positions of the different markers linked to Mt-FTQTL6 including 
MTIC153 and NT6054 that border its confidence interval. 

Additional file 3: Construction of a BAC clone contig map spanning 
Mt-FTQTL6-containing region. Describes the strategy that was followed 
to close the physical gap coinciding with Mt-FTQTL6 and summarizes the 
main results. The text is organized in three sections: (1) State of art; (2) 
Methods and (3) Achievements. 

Additional file 4: Positive PCR-conflrmed BAC clones obtained 
through screening of high-density colony filters with 15 probes. 

Lists the BAC clones from M. truncatula mth2 and Mtf83 libraries that 
were identified by screening high-density colony filters with probes 
corresponding to 15 Mt-FTQTL6-linked markers and were further 
confirmed positive by PCR control. 

Additional file 5: BAC clone contigs partly spanning M. truncatula 
chromosome 6 region harboring Mt-FTQTL6. Depicts the six BAC 
clone contigs that were constructed with BAC clones identified either 
through screening of M truncatula BAC libraries mth2 and Mtf83 
(see Additional file 4) or using an in silico approach. The BAC clone 
contig map is connected to the genetic map of (Vlt-FTQTL6 by means of 
the 15 markers used as screening probes and other additional markers. 
The three BAC clones that were selected for insert sequencing are 
shown. 

Additional file 6: Statistics of 454 BAC sequencing and contig 
assembly. Contains information relative to 454 sequencing of BAC 
clones mth2-92015, mtel-60A22 and mth2-221P20. Statistics on raw, 
cleaned and assembled sequences are given. 

Additional file 7: Collnearlty conservation between Mt-FTQTL6 
Interval and chromosomal regions from sequenced dicotyledonous 
genomes. Illustrates colinearity conservation between Mt-FTQTL6 region 
and different genomic segments from reference dicotyledonous plant 
species including M. truncatula, G. max, Ljaponlcus, P. trichocarpa, 
A. thallana, V. vinlfera and S. lycopersicum. 

Additional file 8: Pairwise comparison of predicted protein 
products of MtCBF2-3;5-10;12-14. Shows amino acid identity and 
similarity between the predicted protein products ot MtCBF2-3:5-IO:l2-l4. 
This information is important to understand the phylogenetic 
relationships between tandemly duplicated CBF/DREBI genes. 

Additional file 9: Amino acid variants within CBF/DREB1 signature 
sequences. Reports the amino acid variants in conserved CBF/DRFBl 
signature motifs PKKP/RAGRxKFxETRHP, DSAWR, A(A/V)xxA(A/V)xxF and 
LWSY that are carried by MtCBF2-3;5-l 1;14. Other CBF/DRFBl proteins 
sharing the same amino acid variant(s) are indicated when possible. 

Additional file 10: Sequencing of MtCBF2-3;5-12;14: primers, PCR 
conditions and sequence analyses. Contains details on the primers 
that were used to isolate and sequence A/JOf2-5;5-/2;/4 from F83005-5 
and DZA045-5. For each of the CBF/DREBI genes, SNP, indel and SSR 
polymorphisms identified between F83005-5 and DZA045-5 together 
with amino acid variations amongst the predicted protein products are 
indicated. 

Additional file 11: Bibliographic review regarding non-CBF/DREBl 
candidate genes for Mt-FTQTL6. Briefly reviews the current knowledge 
on non-Cfif/D/?f6/positional candidate genes for Mt-FTQTL6 and particularly 
their functions. Information is available for MtBAG-l, MtPFRLD, MTR_050s0020, 



MtZFWD, MTR_054s000l and MTR_054sOOI9 whereas no data could be 
obtained for MTR_050sOOI9 and MTR_6g089580. 

Additional file 12: Expression profiles of a subset of M. truncatula 
CBF/DREBI genes according to the Gene Expression Atlas. Shows the 
mean transcript levels from three independent biological replicates of 
MtCBF2, MtCBFS, MtCBFS, MtCBF6, MtCBF/, MtCBF8.1, MtCBF8.2 and MtCBF9 
in four unstressed organs (root, stem, leaf and flower) of the Al 7 
accession. Data were obtained from M. truncatula Gene Expression 
Atlas version 2 [76]. 
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